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Abstract 

Measurements are presented of i?b, the ratio of the bb cross-section to the qq 
cross-section in e + e~ collisions, and the forward-backward asymmetry Ap B at 
twelve energy points in the range ^/s = 130 — 207 GeV. These results are found 
to be consistent with the Standard Model expectations. The measurements are 
used to set limits on new physics scenarios involving contact interactions. 
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1 Introduction 

The ratio i?b = o"(e + e~ — >■ bb)/a(e + e~ — > qq) and Ap B , the forward-backward pro- 
duction asymmetry of bottom quarks in e + e~ collisions, are important parameters in 
precision studies of electroweak theory, and are sensitive probes of new physics. This 
paper presents measurements of and A FB made at centre-of-mass energies (\/s) be- 
tween 130 GeV and 207 GeV. Events containing a bb pair have several characteristic 
features, most notably the presence of secondary vertices, which may be used to select 
a sample enriched in b-decays. A 'b-tag' variable has been constructed for this purpose, 
which exploits the high resolution tracking provided by the DELPHI Silicon Tracker. 
In the asymmetry measurement the hemisphere containing the b-quark has been deter- 
mined using a hemisphere-charge technique. In order to enhance sensitivity to possible 
new physics contributions from high energy scales, all measurements have been made for 
events in which \Js'/s > 0.85, where is the effective centre-of-mass energy after initial 
state radiation. In the Standard Model e + e~ — > bb events are produced by an s-channel 
process propagated by either photon or Z-boson exchange. Over the interval of collision 
energies under investigation the relative strengths of the two contributions evolve so that 
the value of is expected to fall, and that of A FB to rise, slowly with y/s. 

Studies of bb production at collision energies above the Z-pole have been presented by 
other LEP collaborations [1-5]. The results presented here for the energies 130 < yfs < 
172 GeV supersede those of an earlier DELPHI publication [6]. 

Sect. 2 describes the datasets and the aspects of the DELPHI detector relevant for the 
analysis. The event selection is discussed in Sect. 3. The i?b determination is presented in 
Sect. 4 and that of A\ B in Sect. 5. An interpretation of the results within the context of 
both the Standard Model and possible new physics models including contact interactions 
is given in Sect. 6. 

2 Datasets, the DELPHI Detector and Simulation 

LEP 2 operation began in 1995, when around 6 pb _1 of data were delivered at centre- 
of-mass energies of y/s =130 GeV and 136 GeV. In 1996 the collision energy of the 
beams was raised to, and then beyond, the W + W~ production threshold of 161 GeV. 
Each subsequent year saw increasing amounts of integrated luminosity produced at ever 
higher energies, reaching 209 GeV in the year 2000. In total around 680 pb _1 were 
collected by the DELPHI experiment at 12 separate energy points. Note that during 
the 2000 run, operation occurred at a near-continuum of energies between 202 GeV and 
209 GeV. In the present study the data collected during 2000 are divided into two bins, 
above and below 205.5 GeV. Throughout LEP 2 operation collisions were performed 
with unpolarised beams. The mean collision energies for each period of operation and 
the integrated luminosities used in the analysis are summarised in Table 1. More details 
on the LEP collision energy calibration and the DELPHI luminosity determination are 
given in [7] and [8], respectively. 

In addition to the high energy operation, in each year from 1996 onwards LEP also 
delivered 1-4 pb -1 at the Z-pole, in order to provide well understood calibration data for 
the experiments. In this paper the events collected during the calibration running are 
referred to as the 'Z-data', and provide control samples for the high-energy studies. In 
1995 the control sample is taken from the Z-peak data immediately preceeding the switch 
to 130 GeV operation. In 2000 a second set of Z-data was collected in order to provide 
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a dedicated calibration sample for the period in which the DELPHI TPC had impaired 
efficiency (see below). 

A description of the DELPHI detector and its performance can be found in [9,10]. For 
the analyses presented in this paper, the most important sub-detector in DELPHI was 
the Silicon Tracker [11]. The Silicon Tracker was a three-layer vertex detector providing 
measurements in both the views transverse and longitudinal to the beam line, with the 
capabilities to provide effective b-tagging over the polar angle interval of 25° < 9 < 155°, 
where 9 is the angle with respect to the e~ beam direction. End-caps of mini-strip and 
pixel detectors gave tracking coverage down to 9 = 10° (170°). The Silicon Tracker was 
fully installed in 1996 and remained operational until the end of the LEP 2 programme. 
During the 1995 run b-tagging information was provided by the microvertex detector 
described in [12]. 

During the 2000 run, one of the 12 azimuthal sectors of the central tracking chamber, 
the TPC, failed. After the beginning of September 2000 it was not possible to detect the 
tracks left by charged particles in that sector. The data affected correspond to approxi- 
mately one quarter of the total dataset of that year (the 'BTPC period). Nevertheless, 
the redundancy of the tracking system of DELPHI meant that tracks passing through 
the sector could still be reconstructed from signals in the other tracking detectors. A 
modified tracking reconstruction algorithm was used in this sector, which included space 
points reconstructed in the Barrel RICH detector. As a result, the track reconstruction 
efficiency was only slightly reduced in the region covered by the broken sector, but the 
track parameter resolutions were degraded compared with the data taken prior to the 
failure of this sector (the 'GTPC period). 

To determine selection efficiencies and backgrounds in the analysis, events were sim- 
ulated using a variety of generators and the DELPHI Monte Carlo [10]. These events 
were passed through the full data analysis chain. Different software versions were used 
for each year, in order to follow time variations in the detector performance. For the year 
2000, separate GTPC and BTPC sets of simulation were produced. The typical size of 
the simulated samples used in the analysis is two orders of magnitude larger than those 
of the data. 

The e + e~ — ■> ff process was simulated with KK 4.14 [13], interfaced with PYTHIA 
6.156 [14,15] for the description of the hadronisation. For systematic studies, the alter- 
native hadronisation description implemented in ARIADNE 4.08 [17] was used. Four- 
fermion background events were simulated with the generator WPHACT 2.0 [18,19], with 
PYTHIA again used for the hadronisation. 

3 Event Selection 

The analysis was made using charged particles with momentum lying between 0.1 GeV 
and 1.5-(y / i/2), and measurement uncertainty of less than 100%, and having a closest 
approach to the beam-spot of less than 4 cm in the plane perpendicular to the beam 
axis, and less than 4/sin# cm along the beam axis. Neutral showers were used above a 
minimum energy cut, which was 300 MeV for the barrel electromagnetic (HPC) and very 
forward calorimeter (STIC), and 400 MeV for the forward electromagnetic calorimeter 
(FEMC). 

The following requirements were applied to select a pure sample of hadronic events, 
and to ensure that each event lay within the acceptance of the Silicon Tracker: 

• Number of charged particle tracks > 7; 
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• Quadrature sum over each end-cap of energy reconstructed in the forward electro- 
magnetic calorimeter system (STIC + FEMC) < 0.85(^/2); 

• Total transverse energy > 0.2^; 

• Energy of charged particles > 0.1 y/s; 

• Restriction on the polar angle of the thrust axis, 9t, such that | cos#t| < 0.9. 

Data-taking runs were excluded in which the tracking detectors and Silicon Tracker were 
not fully operational. 

In addition to this selection a 'W-veto' was applied to suppress the contamination 
from four-fermion events. The veto procedure consisted of forcing the event into a four- 
jet topology using the LUCLUS [14,15] algorithm and imposing the requirement that 
(-Emin/v/i) ■ a m i n < 4.25°, where E min is the energy of the softest jet, and a m i n the 
smallest opening angle found between all two-jet combinations. This condition is designed 
to distinguish between two-fermion events containing gluon jets, and genuine four-fermion 
background. Less than 40% of four-fermion events survive the hadronic selection and the 
W-veto. 

The analysis is concerned with events produced with an effective centre-of-mass energy 
of the qq system, \fs' \ at or around the collision energy, y/s. The effective centre-of-mass 
energy is reconstructed as in the hadronic analysis reported in [8]. A constrained fit 
is performed, taking as input the observed jet directions as found by the DURHAM 
clustering algorithm [16], imposing energy and momentum conservation, and assuming 
any ISR photon was emitted along the beam line. Radiative returns to the Z are then 
rejected by requiring that the reconstructed value of yj s' / s > 0.85. Contamination from 

events with true values of \J s' / s below this threshold is around 16% at 130.3 GeV and 
reduces to about 6% at 206.6 GeV. 

As a final condition, events with |Qp B | > 1.5 are rejected, where |<5fbI i s one °f the 
event charge variables defined in Sect. 5.1. This selection is applied to exclude badly 
measured events from the asymmetry measurement, and removes around 0.5% of the 
sample. 

The numbers of events passing the high \J s' / s two-fermion hadronic selection at each 
energy point are listed in Table 1, together with the Monte Carlo expectations. The two 
sets of numbers agree well. The background from four-fermion events is estimated to 
be around 9% in the 172.1 GeV dataset, rising to 21% in the 206.6 GeV sample. The 
contamination from t + t~ events is around 0.3%. All other backgrounds are negligible. 

A 'b-tag' variable is used to extract a sub-sample of events enriched in b-quarks from 
the non-radiative qq sample. This variable makes use of three observables, known to 
distinguish between b-quark events and those events with non-b content. In this analysis, 
the three categories of observable considered are: 

• A lifetime variable, constructed from the impact parameters of charged particle 
tracks in each jet; 

• The invariant mass of charged particles forming any secondary vertices that are 
found; 

• The rapidities of charged particles in any secondary vertex, defined with respect to 
the jet direction. 

These properties are used to construct a single event 'b-tag' variable, B tag , of typical value 
between -5 and 10. Events with higher values of this variable are enriched in b-events. 
More information on the b-tagging procedure may be found in [20]. In this analysis a 
cut value of 1 is used for all high energy data sets to select the b-enriched sample; this 
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Table 1: The year of data-taking, mean centre-of-mass energy, integrated luminosity, 
number of events after hadronic selection and W-rejection ('Before b-tag'), and number 
of events after the b-tag. In the year 2000 the numbers in parentheses are those corre- 
sponding to the GTPC sub-sample. Numbers are shown for data and Monte Carlo, where 
for the latter the samples have been scaled to the integrated luminosity of the data and 
Standard Model cross-section values are assumed. 









Before b-tag 


After b-tag 


Year 


[GeV] 


JCdt [pb- 1 ] 


Data 


MC 


Data 


MC 


1995 


130.3 


2.9 


224 


224 


30 


24 




136.3 


2.6 


160 


160 


15 


17 


1996 


161.3 


10.1 


363 


321 


46 


36 




172.1 


10.0 


304 


280 


27 


29 


1997 


182.7 


53.1 


1351 1284 


117 137 


1998 


188.6 


156.8 


3567 


3541 


365 


379 


1999 


191.6 


25.8 


563 


565 


68 


57 




195.5 


76.2 


1629 


1597 


164 


159 




199.5 


83.0 


1651 


1670 


184 


162 




201.7 


40.6 


807 


799 


88 


77 


2000 


204.8 (204.8) 


82.8 (76.1) 


1538 (1411) 1572 (1447) 


144 (131) 147 (137) 




206.6 (206.6) 


136.4 (84.7) 


2510 (1586) 2536 (1581) 


240 (167) 233 (148) 




Total 


680.3 


14667 


14549 


1488 


1457 



selection has a typical efficiency for bb events of around 65%, but only 2.5% for cc events 
and 0.3% for light quark events. The numbers of events passing the b-tag are listed in 
Table 1. Here the Monte Carlo numbers do not include the correction factors discussed 
in Sect. 4. 



4 Measurement of i? D 

4.1 Procedure and Calibration with Z Data 

For each energy point i? D is determined through the following relation: 
_/V D - N 4i 

— g — _ 4f = RhCb^b + R c c c e c + e uds (1 — c c R c — c h R b ). (1) 

-''total -''total 

Here -^otai (tag) an< ^ -^totai(tag) are the number of events in the data, and the estimated 
four-fermion background respectively, before (after) the application of the b-tag cut; R c is 
directly analogous to i? D , but defined for cc events; and e&, e c and e uds are the efficiencies 
of the b-tag cut applied to b, c and light quark events respectively, q, and c c are correction 
factors, which account for the fact that the effective values of _R D and R c are modified by 
the hadronic selection, and that there is some contamination from initial state radiative 
production in the sample, the fraction of which can in principle be different for each quark 
type, and therefore changes with the application of the b-tag. Simulation indicated that 
these correction factors lie within 1-2% of unity. 

The efficiency and expected background were determined primarily from Monte Carlo, 
and cross-checked, where possible, from the data themselves. Figure 1 shows the distri- 
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bution of the b-tag variable, B tag , in data and simulation for each dataset. In these plots 
the 2000 data have been divided between GTPC and BTPC operation, and the 1995 
and 1996 data have been combined. In general, reasonable agreement can be seen for all 
years in the region around and above the cut position of B tag = 1.0, with worse agreement 
for the background-dominated region below the cut. (The implications of this imperfect 
background description are assessed below.) 

The running at the Z-pole in each year provides a control sample which may be used 
to calibrate the simulation. The value of i?b at the Z-pole is well known from LEP 1 [21]. 
This value has been compared with the results obtained from applying expression (1) to 
each sample of Z-calibration data. Figure 2 shows the distribution of B tSug for Z-calibration 
data of the 2000 GTPC period, together with that of the corresponding simulation. The 
b-tag variable has a mild dependence on the collision energy. In order to make the Z-data 
study as relevant as possible to the high energy measurements, the cut value was placed 
at Stag = 0.6 for these data, which gives a similar efficiency to the value used at high 
energy. The analysis returned a value of i?b which was similar for all datasets apart from 
1998, with a mean that was (4.1 ± 1.2)% higher in relative terms than the world average 
result. The value found for 1998 was (4.2 ± 1.4)% lower than the world average. 

The offset in the measurement of R h with the Z-data can be caused by imperfections 
in simulating the response of the detector to the b events, the background or to both. 
(Effects arising from uncertainties in the knowledge of the B and D decay modelling 
have been accounted for and found to be small.) In order to distinguish between these 
possibilities, a fit was performed to the B tag distribution of the Z-data in the background 
enriched region around the cut value (0 < B tag < 2.5), taking the shapes of the signal 
and background from the simulation and fitting their relative contributions. The results 
returned background scaling factors with respect to the simulation which varied between 
around 0.9 and 1.2, depending on the year, with a relative precision of better than 5%. 
After allowing for these corrections, the remaining, and most significant, cause for the 
offset was attributed to an incorrect estimate of the b-tagging efficiency. 

A fit was performed to the background level in the high energy data, identical to 
that made with the Z-running samples. Compatible results were obtained within ±10%. 
For the high energy R h extraction, therefore, these Z-pole determined scaling factors were 
applied to the cc and uds background, with this 10% uncertainty assigned as a systematic 
error, uncorrelated between years. The same factors were applied to the four-fermion 
background, but with twice the systematic uncertainty, as this background component is 
not present in the Z-data. Finally, the b-tagging efficiency was corrected by the amount 
indicated from the low energy study, with half of this correction taken as an uncertainty, 
to account for any variation with energy. The correction factor varied between 0.959 in 
1998 and 1.045 for the highest energy point of 2000. Given the very similar nature of the 
offset seen in the Z-pole study for all years apart from 1998, the uncertainty was taken 
as correlated for these datasets. 

The calibration procedure was repeated under different conditions and assumptions, 
for example using the same B tSug cut value for Z-pole and high energy data, and using 
an absolute offset rather than a factor to correct the efficiency. In all cases compatible 
results were obtained. 

Table 2 shows the post b-tag sample composition at each energy point, after applying 
the various corrections factors and assuming the Standard Model production fractions. 
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Figure 1: The variable B tag plotted for all datasets. The standard analysis has a cut at 
-Stag = 1- The insets show a zoom of the b-enhanced region on a linear scale. 



7 




Figure 2: The variable B tag for the 2000 GTPC Z-data. The inset shows a zoom of the 
b-enhanced region on a linear scale. 



Table 2: The percentages of each event category making up the sample after the cuts on re- 
constructed \J s'/ s and B tag , for each energy. In the case of bb events the division between 

high and low true \J s' / s is indicated. For the non-bb final states the low yfs' component 
is included in the category definition. (Note that for the energies y/s =130.3-172.1 GeV, 
the cc and uds background contributions have uncertainties of around 0.5%, due to finite 
Monte Carlo statistics.) 



yfi [GeV] 


bb rrr^n., 

^/s'/s>0.85 


bb nr nBK 


cc 


uds 


4-fermion 


130.3 


79.7 


15.8 


3.7 


0.7 





136.3 


77.8 


17.9 


2.9 


1.4 





161.3 


83.9 


10.6 


4.4 


0.5 


0.6 


172.1 


82.3 


8.4 


4.9 


1.6 


2.7 


182.7 


82.1 


7.7 


5.1 


1.5 


3.6 


188.6 


81.8 


7.0 


5.6 


1.5 


4.2 


191.6 


83.1 


6.9 


4.9 


1.2 


4.0 


195.5 


82.7 


6.7 


4.9 


1.5 


4.2 


199.5 


82.9 


6.5 


4.9 


1.3 


4.4 


201.7 


82.6 


6.4 


4.8 


1.5 


4.6 


204.8 


81.7 


6.2 


5.3 


1.5 


5.2 


206.6 


82.1 


6.0 


5.0 


1.6 


5.3 
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4.2 Systematic Uncertainties in Modelling of Physics 
Processes 

The stability of the results was studied with respect to uncertainties in the knowledge 
of important properties of B and D production and decay, and other event characteristics 
relevant to the b-tag. The variation in the parameter values was implemented by re- 
weighting Monte Carlo events to the modified distribution. 

• b and c fragmentation: Simulated bb and cc events at high energy had their 
Peterson fragmentation parameters [22] varied in the range corresponding to the 
uncertainties in the mean scaled energy of weakly decaying b and c hadrons in Z 
decays [21]. 

• b and c decay multiplicity: The charged b decay multiplicity was allowed to 
vary in the range 4.955 ± 0.062 [21] and that of D mesons was varied according 
to [21,23], with a ±0.5 uncertainty assigned to the charged multiplicity of c baryon 
decays. 

• b and c hadron composition: The proportions of weakly decaying b and c hadrons 
were varied according to the results reported in [24] and [25] respectively. 

• b and c hadron lifetime: The b and c hadron lifetimes were varied within their 
measured range [24]. In the b hadron case this was 1.576 ± 0.016 ps. 

• gluon splitting to heavy quarks: The rate of gluon splitting to bb and cc per 
hadronic event was varied in the range (0.254 ± 0.051)% and (2.96 ± 0.38)% respec- 
tively [21]. 

• Kg and A production: The rate of Kg and A hadrons was varied by ±5%, 
consistent with [26,27]. 

For each property in turn, the value of Rb was recalculated using the re-weighted simu- 
lation as input and the observed change taken as the systematic uncertainty. The results 
for the 188.6 GeV and 206.6 GeV energy points are shown in Table 3, with the total un- 
certainty corresponding to the sum in quadrature of the individual components. Similar 
behaviour was observed for the other energy points. 

4.3 Summary of Systematics and Results 

The relative systematic uncertainties on R h are summarised in Table 4. In addition 
to those components already discussed, contributions are included which arise from the 
finite size of the Monte Carlo simulation sample, and from the effect of the uncertainty 
in the residual radiative contamination in the analysis. Studies on the resolution of the 



ys'/s reconstruction indicated that this background was understood to the level of 10%. 
It can be seen that the dominant source of systematic uncertainty is that coming from 
the comparison with the Z-data. 

The results for Rb are given in Table 5, together with the statistical and systematic 
uncertainties. The correlation matrix for these results can be found in Appendix A. For 
each of the two energy points of the year 2000 the results for the GTPC and BTPC period 
are found to be compatible and are thus combined into a single value. No variation of R c 
is considered in the systematic uncertainty, but the dependence of Rb on this quantity, 
AR h /(R c - R^ M ), is tabulated explicitly. 

The internal consistency of the measured Rb results may be studied, under the as- 
sumption that any dependence of the true value on collision energy can be neglected. 
The pull distribution of (Rb~ < Rb >)/er is found to have a spread of 1.2, with the most 
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Table 3: Fractional systematic uncertainties on R b associated with physics modelling for 
two illustrative energy points. Values are given in percent. 





Energy point 


Uncertainty bource 


188.6 GeV 206.6 GeV 


b fragmentation 


0.2 


0.2 


b decay multiplicity 


0.5 


0.7 


b hadron composition 


0.2 


0.2 


b hadron lifetime 


0.2 


0.3 


c fragmentation 


0.1 


0.1 


c decay multiplicity 


0.3 


0.2 


c hadron composition 


0.2 


0.2 


c lifetime 


0.1 


0.1 


g — > bb 


0.1 


0.1 


g -> cc 


0.1 


< 0.1 


Kg and A production 


0.2 


0.3 


Total 


0.8 


0.9 



Table 4: The fractional systematic uncertainty, in percent, on R h , energy point by energy 
point. 



[GeV] 


Z Comparison 


Modelling 


4-fermion MC Stats 


Rad. Bckgd. 


Total 


130.3 


1.7 


1.1 


/ 


2.4 


0.5 


3.2 


136.3 


1.8 


1.1 


/ 


2.9 


0.4 


3.6 


161.3 


1.6 


1.1 


0.1 


1.9 


0.1 


2.7 


172.1 


1.8 


1.1 


0.5 


2.2 


0.1 


3.1 


182.7 


2.1 


1.0 


0.9 


0.4 


0.1 


2.5 


188.6 


2.1 


0.8 


0.8 


0.4 


0.1 


2.4 


191.6 


2.2 


0.8 


0.8 


0.5 


0.3 


2.5 


195.5 


2.2 


0.9 


0.9 


0.5 


0.1 


2.6 


199.5 


2.2 


1.0 


0.9 


0.5 


0.2 


2.6 


201.7 


2.2 


0.9 


0.9 


0.4 


0.2 


2.6 


204.8 


2.1 


0.9 


1.2 


0.4 


0.1 


2.6 


206.6 


2.4 


0.9 


1.1 


0.3 


0.2 


2.8 
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Table 5: The results for at each energy point. Also given are the dependences of i?b on 
R c , and the values for the latter fraction assumed in the analysis [28]. For convenience, 
the corresponding Standard Model expectations for R h are included. 



,/s [Op VI 


Rb 




C"stat 




IT 

'-'syst 




pSM 
n c 


r< 


SM 










b 


1 on o 
130.3 





228 


i 

± 





041 


i 

± 


0.00 ( 


-0 


027 


O.zzO 





lo6 


136.3 





153 


± 





041 


± 


0.006 


-0 


023 


0.226 


o 


182 


161.3 





183 


± 





029 


± 


0.005 


-0 


023 


0.244 





170 


172.1 





127 


± 





028 


± 


0.004 


-0 


023 


0.249 





167 


182.7 





127 


± 





013 


± 


0.003 


-0 


032 


0.253 





165 


188.6 





166 


± 





009 


± 


0.004 


-0 


035 


0.255 





164 


191.6 





194 


± 





024 


± 


0.005 


-0 


032 


0.256 





163 


195.5 





161 


± 





013 


± 


0.004 


-0 


031 


0.258 





163 


199.5 





187 


± 





014 


± 


0.005 


-0 


031 


0.258 





162 


201.7 





183 


± 





020 


± 


0.005 


-0 


030 


0.259 





162 


204.8 





156 


± 





014 


± 


0.004 


-0 


031 


0.259 





161 


206.6 





163 


± 





011 


± 


0.005 


-0 


029 


0.260 





161 



outlying entry arising from the measurement at y/s = 183 GeV, which is 2.7 a below the 
mean. 

The stability of the results has been examined when changing the value of the b-tag 
cut. The cut position was tightened to a value of B tag =2.5 in the high energy data, and 
-B ta g=2.1 in the Z-data, and i?b re-evaluated at each energy point. Under this selection 
the event samples halve in size, but the non-bb background is reduced by almost a factor 
of three. No statistically significant change in result was observed with respect to the 
standard selection for any energy point in isolation, nor for all energy points averaged 
together, indicating that the background levels and efficiency are well understood for 
both selections. 

The results for i?b are compared with the Standard Model expectations and interpreted 
in the context of possible new physics contributions in Sect. 6. 



5 Measurement of Ap B 
5.1 Procedure 

For the non-radiative bb events selected in this study, the expected form of the differ- 
ential cross-section is given by: 

, dC7h n oc 1 + cos 2 # b + !,4k B cos0b, (2) 
a cos t^b 3 

where 9^ is the polar angle the b-quark makes with the initial e~ direction. 

The analysis presented in this paper is based on an unbinned likelihood fit to expres- 
sion (2), and hence requires knowledge of 6 h cc , which is the event-by-event value of 9 h as 
reconstructed in DELPHI. This reconstruction is performed using the thrust axis and a 
hemisphere charge technique. Each event is divided into two hemispheres by the plane 
perpendicular to the thrust axis that contains the nominal interaction point. Simulation 
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shows that for non-radiative events the thrust axis is a good approximation to the di- 
rection of emission of the initial bb pair. Then the 'hemisphere charges' Qp and Q B are 
calculated for the forward and backward hemispheres. Qp is defined: 

Qr = E,|P,-T|"' (3) 

where p« and g« are the momentum and charge of particle i, T is the thrust 
is an empirical parameter, and the sum runs over all charged particle tracks for which 
Pi ■ T > 0. Qb is defined in an analogous manner with the requirement that • T < 0. 
The information from both hemispheres may be combined into two event variables: 

Q± B = Q F ±Q B . (4) 

The sign of Qp B is sensitive to whether the b-quark was emitted in the forward or back- 
ward hemisphere. The value of k in equation (3) is tuned to maximise this discrimination, 
and is set to 0.5. Figure 3 (a) shows Qp B , plotted for all data. There is a small, but 
significant negative offset, indicating that the b-quark is preferentially emitted in the for- 
ward hemisphere. Qp B has no sensitivity to the initial b-quark direction, but provides a 
quantity which can be compared between data and simulation, with a width that reflects 
the resolution of the method. Qp B is plotted in Fig. 3 (b), together with the correspond- 
ing quantity from the simulation. As expected, it is centred on zero. The distribution is 
marginally wider in data than in the Monte Carlo. 




Figure 3: Distribution of the two event charge variables for all data after b-tag cut. (a) 
shows the charge asymmetry between the two hemispheres, Qp B - (b) shows the sum of 
the hemisphere charges, Qp B . Also shown are the expectations from the simulation, which 
are generated with the Standard Model values for the asymmetries of each component. 



The cosine of the reconstructed b-quark direction is then given by: 

X ee cos6£ oc = - sign(Qp B ) • | cos# T |, (5) 

where Qt is the polar angle of the thrust axis. The distribution of cos0£ ec is shown 
in Fig. 4 (a), for the full LEP 2 dataset, plotted for events where \Qrp B \ > 0.1. The 
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asymmetry which is observed is an underestimate of the real asymmetry, both because of 
'mistags' and because of background contamination. Detector inefficiencies also distort 
the distributions, particularly in the forward and backward regions. Mistags are events 
in which the sign of Qp B does not give the correct b-quark direction. Mistags dilute the 
true asymmetry by a factor D — (1 — 2u), where u is the probability of mistag. Note 
that uj has a dependence on the absolute value of Qp B - For example, simulation indicates 
that for the ensemble of high energy data the mistag rate has a value of uj = 0.45 for 
events where \Qp B \ < 0.1, and uj = 0.27 in the case when \Qp B \ > 0.1, falling to uj = 0.17 
when \Qyb\ > 0.36. Figure 4 (b) shows the same data after correction for background 
contamination, detector inefficiency and mistags, and the corresponding distribution for 
the Z-data. It is apparent that the high energy data exhibit an asymmetry significantly 
higher than that of the Z-data, which have a value consistent with that measured at 
LEP 1 [21]. 
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Figure 4: The observed angular distribution for all data after b-tag cut and the require- 
ment |QfbI > 0-1- ( a ) shows the raw distribution of events with respect to cosd™ c to- 
gether with the expectations from simulation, generated with the Standard Model values 
for the asymmetries of each component, (b) shows the differential cross-section (nor- 
malised to the total cross-section within the acceptance) with respect to cos6'£ or , where 
#£ or is the b-quark direction after correction for wrong flavour tags, non-uniform accep- 
tance efficiency and background. Also shown is the corresponding distribution for the 
LEP 2 Z-data. The superimposed curves are fits to the form of the expected differential 
cross-section. 



Optimal sensitivity to Ap B is achieved through performing a maxmimum likelihood 
fit, taking as the probability density function the expected differential cross-section of 
equation (2). At each energy point, the measured asymmetry A™ B as is determined by 
maximising the following expression: 

ln£ = J> (l + (A,) 2 + |aS§-X,) , (6) 
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where the sum runs over all events. Mistags and contamination are accounted for by 
writing 

Y.f,l\r\r (7) 

j 

Here the sum runs over the five categories of event type in the sample: signal, radiative bb 
contamination, cc, light quark and four-fermion. Each category enters with a proportion 
fj, as given by the values in Table 2, with a true asymmetry Aj and dilution factor Dj, 
where Aj for the signal category is equivalent to Ap B . For the purposes of accounting for 
the background in the fit, equation (2) is an adequate description of the distribution of 
radiative and four-fermion events. The dilution factors are determined from simulation, 
and the asymmetries of the background processes are set to their Standard Model ex- 
pectations. In order to exploit the dependence of the mistag probability on the absolute 
value of the charge asymmetry, all events are used, but the dilutions and event fractions 
are evaluated in four bins of |<5fbI an d included in the fit accordingly. 

The fit procedure has been tested on a large ensemble of simulated experiments, and 
found to give unbiased results with correctly estimated uncertainties. It has also been 
applied to the Z-data. Averaged over all datasets, the measured asymmetry minus that 
value determined at LEP 1 [21] is found to be -0.01 ± 0.01. 

5.2 Results and Systematic Uncertainties 

The most important source of systematic uncertainty in the asymmetry measurement 
is associated with the knowledge of the performance of the charge asymmetry variable. 
There are three significant contributions to this uncertainty: 

• Detector Response: The distribution of track multiplicity as a function of mo- 
mentum has small differences between data and Monte Carlo both at high and low 
momentum, which may be attributed to an imperfect modelling of the track recon- 
struction in the simulation. Tracks were re-weighted in the simulation in order to 
establish the effect on the mistag rate. Similar studies were conducted to under- 
stand the consequences of differences in the momentum resolution between data and 
Monte Carlo. Finally, the width of the Qp B distribution was artificially increased in 
the simulation, to match that of the data, by adjusting the value of the n parameter 
in the analysis of the simulation alone, and the effect on Qp B was determined. 

• Hadronisation: An alternative Monte Carlo data set of events based on ARI- 
ADNE [17] was used to assess the robustness of the estimation of the mistag rate 
with respect to the description of the hadronisation process used in the simulation. 

• Monte Carlo Statistics: The limited amount of simulation data available intro- 
duces a non-negligible statistical uncertainty in the knowledge of the mistag rate. 

Additional possible sources of measurement bias related to the mistag have been consid- 
ered, for example whether any significant angular dependence exists in the value of the 
dilution. These effects were found to have negligible impact on the results. 

In addition to these studies, systematic uncertainties were evaluated arising from the 
same three sources that were considered in the measurement, namely the uncertainty 
associated with the sample composition as assessed from the Z-data; the uncertainty 
in the level of the 4-fermion background; and the uncertainty in the modelling of the 
physics processes (apart from hadronisation). The modelling systematic here includes a 
component arising from the uncertainty in the knowledge of the b-mixing parameter x- 
This was varied within the range 0.128 ± 0.008, following the evaluation reported in [24]. 
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Table 6: Systematic uncertainties on A^ B for two illustrative energy points. 





Energy point 


Uncertainty Source 


188.6 GeV 206.6 GeV 


Detector Response 


0.054 


0.038 


Hadronisation 


0.027 


0.025 


MC Statistics 


0.016 


0.011 


Z Comparison 


0.008 


0.004 


Modelling 


0.008 


0.008 


QCD Correction 


0.018 


0.018 


4-fermion 


0.003 


0.006 


Radiative background 


0.004 


0.004 


Total 


0.066 


0.051 



Table 7: The results for Ap B at each energy point, together with the Standard Model 
expectation [28]. 



x/i [GeV] 


^FB 




°stat 




^"syst 


A 


b,SM 
FB 


130.3 


0.569 


± 





507 


± 





112 





473 


136.3 


0.447 


± 





615 


± 





117 





496 


161.3 


1.344 


± 





346 


± 





097 





550 


172.1 


0.407 


± 





523 


± 





099 





564 


182.7 


-0.120 


± 





245 


± 





102 





575 


188.6 


0.703 


± 





157 


± 





066 





579 


191.6 


0.391 


± 





304 


± 





049 





582 


195.5 


0.875 


± 





221 


± 





060 





584 


199.5 


0.602 


± 





185 


± 





052 





587 


201.7 


0.756 


± 





298 


± 





055 





588 


204.8 


0.718 


± 





252 


± 





061 





590 


206.6 


0.108 


± 





180 


± 





051 





591 



A further uncertainty is assigned to account for the fact that QCD corrections to the final 
state, in particular gluon radiation, modify the asymmetry. The size of this effect has 
been estimated using ZFITTER [28] to be 0.018. In practice the selection cuts disfavour 
events with hard gluon radiation and thus will suppress this correction. In this study, 
however, the full effect is taken as an uncertainty, fully correlated between energy points. 
Finally, a systematic error is added to account for the uncertainty in the knowledge of 
the residual radiative bb contamination in the sample. 

Table 6 lists the systematic uncertainties for the 188.6 GeV and 206.6 GeV energy 
points. The total is the sum in quadrature of the uncorrelated component uncertainties. 
The results for A^ B , including statistical and systematic uncertainties, are shown in 
Table 7. The correlation matrix for these results can be found in Appendix A. Both 
the statistical uncertainty and certain components of the systematic uncertainty have a 
dependence on the absolute value of the asymmetry. The uncertainties shown have been 
evaluated assuming the Standard Model value. 
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The self-consistency of the results may be assessed assuming that any dependence of 
the true value of Ap B on the collision energy can be neglected. The pull distribution of 
(Ap B — < Ap B >)/cr is found to have a spread of 1.5. The outliers contributing to this 
larger than expected width are the dataset at 161.3 GeV, which has an asymmetry which 
is 2.3 a higher than the mean, and the samples at 182.7 GeV and 206.6 GeV, which 
have asymmetries that are low by 2.7 and 2.4 a respectively. The 206.6 GeV dataset is 
made up of events accumulated during both the GTPC and BTPC running; the values of 
the asymmetry and associated statistical uncertainties are found to be 0.087 ± 0.218 and 
0.152 ± 0.318, and hence consistent, for the two periods. All asymmetries have been re- 
evaluated with a more severe b-tag cut of 2.5, as was done for the R\> analysis. Averaged 
over all data points the asymmetry is found to shift by —0.008 ± 0.052 with respect 
to the central values reported in Table 7. The shifts for the 161.3 GeV, 182.7 GeV and 
206.6 GeV datasets are 0.019±0.209, -0.278±0.191 and -0.043±0.162 respectively. The 
magnitudes and signs of these changes do not suggest that there is any significant problem 
with the understanding of the background level and behaviour. Further cross-checks were 
performed in which the fit was restricted to high values of |<5fbI an d where alternative 
methods, such as a binned least-squared fit, were used to determine the asymmetry. 
Again, no significant changes were observed in the results, in particular those of the three 
outlying points. 



6 Interpretation 

The results for R h from Sect. 4.3 and those for Ap B from Sect. 5.2 have been compared 
against the Standard Model expectations, as calculated by ZFITTER [28] with final state 
radiation effects included. The measurements and the expectations are shown in Figs. 5 
and 6, for i?b and Ap B respectively. The mean values of the differences between the 
measurements and the Standard Model expectations have been evaluated using both the 
statistical and systematic uncertainties, and taking full account of all correlations. The 
results of this computation are presented in Table 8. In both cases it can be seen that 
the measurements agree reasonably well with the Standard Model. When all data points 
are combined, the relative precision of the i?b measurements is 3.3% and the overall 
uncertainty on the Ap B measurements is 0.083. These results are the most precise yet 
obtained for the two parameters at LEP 2 energies. 

Contact interactions between initial and final state fermionic currents provide a rather 
general description of the low energy behaviour of any new physics process with a charac- 
teristic energy scale. The results of the i?b and Ap B analyses have been compared with a 
variety of contact interaction models. Following reference [29] the contact interactions are 
parameterised in the same manner as explained in [8], in which an effective Lagrangian 
of the form: 

C eff = | E ViAl^ibj^bj, (8) 

i,j=L,R 

is added to the Standard Model Lagrangian. Here g 2 /^ is taken to be 1 by convention, 
rjij = ±1 or 0, A is the energy scale of the contact interactions, and (bj) are left or 
right-handed electron (b-quark) spinors. By assuming different helicity couplings between 
the initial-state and final-state currents and either constructive or destructive interference 
with the Standard Model (according to the choice of each rjij) a set of different models 
can be defined from this Lagrangian [30]. The values of rjij for the models investigated 
in this study are given in Table 9. 



1G 




Figure 5: The measured values (points) of R h and the Standard Model predictions 
(curve) [28] plotted against y/s. The error bars give the total measurement uncertainties. 




Figure 6: The measured values (points) of Ap B and the Standard Model prediction 
(curve) [28] plotted against y/s. The error bars give the total measurement uncertainties. 
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Table 8: Results of the fit for the mean value of the difference between the measured values 
and the Standard Model predictions [28], for both i?b and Ap B . The first uncertainty is 
statistical, and the second uncertainty is systematic. 



Measurement 


< (Meas - SM) > 


< > [GeV] 


X 2 /ndf (Prob.) 


R h 

^FB 


-0.0016 ±0.0044 ±0.0031 
-0.091 ±0.072 ±0.041 


191.9 
192.2 


17.9/11 (8%) 
20.8/11 (4%) 



In fitting for the presence of contact interactions a new parameter e = 1/A 2 is defined, 
with e = being the limit that there are no new physics contributions. The region e > 
represents physical values of 1/A 2 in models in which there is constructive interference 
with the Standard Model, while the region e < represents physical values for the 
equivalent model with destructive interference. Least squared fits have been made for 
the value of e assuming contact interactions from each model listed in Table 9. All R h and 
Ap B data have been used, taking account of the correlations between the measurements. 
In this fit, the R\, results have been re-expressed as absolute cross-sections, making use 
of the qq cross-section results found in [8]. 

The results of the contact interaction fits are shown in Table 10. The data show no 
evidence for a non-zero value of e in any model, and the table lists the 68% allowed 
confidence level range for the fits to this parameter. Also shown are the corresponding 
95% confidence level lower limits for the contact interaction scale, allowing for positive 
(A + ) and negative (A~) interference with the Standard Model. These limits are in the 
range 2-13 TeV, with the most stringent for the VV, AA and V0 models. 

Table 9: Choices of rjij for different contact interaction models. 



Model 


Vll 


Vrr 


Vlr 


Vrl 


LL ± 


± 1 











RR ± 





± 1 








VV ± 


± 1 


± 1 


± 1 


± 1 


AA ± 


± 1 


± 1 


T 1 


T 1 


LR ± 








± 1 





RL ± 











± 1 


V0 ± 


± 1 


± 1 








A0 ± 








± 1 


± 1 



7 Conclusions 

Analyses of the ratio of the bb cross-section to the hadronic cross-section, i?b, and 
the bb forward-backward asymmetry, Ap B , have been presented for non-radiative pro- 
duction, defined as \J s' / s > 0.85, at 12 energy points ranging from y/s = 130.3 GeV 
to y/s = 206.6 GeV. The relative uncertainties of all i?b measurements is 3.3%, and the 
uncertainty on the mean value of A^ B for all measurements is 0.083, making these results 
the most precise yet obtained for the two parameters at LEP 2 energies. The results are 
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Table 10: Limits of contact interactions coupling to bb. The 68% C.L. range is given for 
e, while 95% C.L. lower limits are given for A^. 



Model 


e (TeV- 2 ) 


A" (TeV) 


A+ (TeV) 


LL 


[-0.0019, 0.0097] 


10.2 


8.4 


RR 


[-0.1947, 0.0172] 


2.2 


5.7 


VV 


[-0.0021, 0.0076] 


10.6 


9.5 


AA 


[-0.0012, 0.0060] 


12.9 


10.7 


LR 


[-0.1029, 0.0234] 


2.9 


4.7 


RL 


[-0.0161, 0.1687] 


5.8 


2.4 


VO 


[-0.0014, 0.0069] 


12.0 


9.9 


AO 


[-0.0163, 0.0630] 


5.3 


3.7 



found to be compatible with those of other experiments [1-5] and are consistent with 
Standard Model expectations. Limits have been derived on the scales of contact interac- 
tions, and are found to lie in the range 2-13 TeV, depending on the chirality structure of 
the new physics contribution. 
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A Correlation Matrices 

The correlation matrices for the i?b and Ap B results are given in Tables 11 and 12 
respectively. The correlations between R h and Ap B are negligible. 



Table 11: Correlation matrix for R h results. 



[GeV] 


130 


136 


161 


172 


183 


189 


192 


196 


200 


202 


205 


207 


130 


1.00 


0.01 


0.01 


0.01 


0.02 


0.01 


0.02 


0.03 


0.03 


0.02 


0.03 


0.04 


136 




1.00 


0.01 


0.01 


0.02 


0.01 


0.01 


0.02 


0.02 


0.01 


0.02 


0.02 


161 






1.00 


0.01 


0.02 


0.01 


0.02 


0.03 


0.03 


0.02 


0.03 


0.04 


172 








1.00 


0.02 


0.01 


0.01 


0.02 


0.03 


0.02 


0.02 


0.03 


183 










1.00 


0.01 


0.04 


0.05 


0.06 


0.04 


0.05 


0.07 


189 












1.00 


0.01 


0.01 


0.02 


0.01 


0.01 


0.02 


192 














1.00 


0.06 


0.06 


0.04 


0.04 


0.06 


196 
















1.00 


0.10 


0.07 


0.06 


0.09 


200 


















1.00 


0.07 


0.07 


0.10 


202 




















1.00 


0.05 


0.07 


205 






















1.00 


0.11 


207 
























1.00 



Table 12: Correlation matrix for results. 



[GeV] 


130 


136 


161 


172 


183 


189 


192 


196 


200 


202 


205 


207 


130 


1.00 


0.04 


0.04 


0.02 


-0.06 


0.07 


0.02 


0.04 


0.04 


0.03 


0.04 


0.05 


136 




1.00 


0.03 


0.02 


-0.04 


0.05 


0.02 


0.03 


0.03 


0.02 


0.03 


0.03 


161 






1.00 


0.05 


-0.06 


0.07 


0.03 


0.05 


0.04 


0.03 


0.04 


0.05 


172 








1.00 


-0.04 


0.04 


0.02 


0.03 


0.03 


0.02 


0.03 


0.03 


183 










1.00 


-0.12 


-0.04 


-0.08 


-0.06 


-0.05 


-0.07 


-0.08 


189 












1.00 


0.05 


0.09 


0.08 


0.06 


0.08 


0.10 


192 














1.00 


0.04 


0.04 


0.03 


0.03 


0.04 


196 
















1.00 


0.06 


0.05 


0.05 


0.06 


200 


















1.00 


0.05 


0.05 


0.06 


202 




















1.00 


0.04 


0.04 


205 






















1.00 


0.06 


207 
























1.00 
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